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Abstract 

A quantitative analysis of = \{g — 2)^ within the framework of Supergravity 
Grand Unification and radiative breaking of the electro-weak symmetry is given. 
It is found that a* USY is dominated by the chiral interference term from the light 
chargino exchange, and that this term carries a signature which correlates strongly 
with the sign of /j,. Thus as a rule a^ USY > for /x > and a^ USY < for /j < with 
very few exceptions when tan/3 ~ 1. At the quantitative level it is shown that if the 
E821 BNL experiment can reach the expected sensitivity of 4 x 10 -10 and there is a 
reduction in the hadronic error by a factor of four or more, then the experiment will 
explore a majority of the parameter space in m,o — rag plane in the region mo ^ 400 
GeV, m g ^ 700 GeV for tan/3 J> 10 assuming the experiment will not discard the 
Permanent address 
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Standard Model result within its 2a uncertainty limit. For smaller tan/3, the SUSY 
reach of E821 will still be considerable. Further, if no effect within 2a limit of the 
Standard Model value is seen, then large tan/? scenarios will be severely constrained 
within the current naturalness criterion, ie., mo,mg ^ 1 TeV. 



1 Introduction 

The high level of experimental accuracy of the measured value of the anomalous magnetic 
moment of the muon [ij has provided verification to several orders in the perturbation 
expansion of quantum electrodynamics(QED) as well as put constraints on physics 
beyond the Standard Model (SM) ||. Further the E821 experiment currently underway 
at Brookhaven is expected to improve the accuracy over the previous measurement by 
a factor of 20 0]. Simultaneously, it is expected that improved analyses of existing 
data on (e + e~ — > hadrons) as well as new data from ongoing experiments in VEPP- 
2M || together with future experiments in DA$NE |J, BEPC etc. will reduce the 
uncertainty in the hadronic contributions to a significant level so as to allow for a test 
for the first time of the elctro-weak corrections in the Standard Model. It was pointed 
out in Refs || |] ( see also Refs |10| for a discussion of previous work and Ref for 



a discussion of recent work ) that supersymmetric corrections to (g — 2)^ are in general 



the same size as the electro- weak corrections in supergravity grand unification [12]. Thus 
improved experiments designed to test the Standard Model electro-weak corrections can 
also provide a probe of supersymmetric contributions. 



2 



Table 1: Contributions to x 10 



Nature of Contribution 


P A QTT A 

priori/ -tV 


PACT? R 


Q.E.D. to 0(a/n) 5 




11658470.8(0.5) 


11658470.8(0.5) 


Kinoshita et al. [ 




i 






Hadronic vac. polarization 


1C\K 't( r 7 ^P.\ 
( \JD.Z\ (.00) 


mo qk/i c op.\ 


to 0(a/7i) 2 






Martinovic & Dubnicka 18 


Eidelman & Jegerlehner 19, 20 


Hadronic vac. polarization 


-9.0(0.5) 


-9.0(0.5) 


to 0(a/n) 3 










Kinoshita & Marciano || 






Light by light hadronic 


U.o^U.y ) 


U.o^U.y J 


amplitude 










Bijnens et al. [21 










Total hadronic 






097.0(7.0) 


694.2(15.3) 


Electro-weak one 


-loop 


ly.o 


ly.o 


Fujikawa et al. [23 








Electro-weak 2-fermion loops 


— z.oyv.o ) 




Czarnecki et al. 
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Electro-weak 2-b 


oson loops 


(-2.0 + 0.045 x R b ) 


(-2.0 + 0.045 x R b ) 


Czarnecki et al. 
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Total electro-weak 




(15.2(0.3) + 0.045 x R b ) 


(15.2(0.3) + 0.045 x R b ) 


upto 2-loops 










Total (with R b = 


0) 


11659183(7.6) 


11659180(15.3) 
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The analyses of Refs [0, ||, however, were done without using the constraints of radiative 



breaking of the electro-weak symmetry |L3[ and without the benefit of the recent high 



precision LEP data |14|, 15] to constrain the coupling constants ||16|| . The purpose of the 



present work is to include these features in the analysis. Additionally, we investigate the 
effects on the results due to b — > s'j constraint |T7[, and dark matter constraint. We shall 



find that the expected accuracy of a /t = ~(g — 2)^, in E821 BNL experiment would either 
allow for supergravity grand unification effects to be visible in the BNL experiment, or 
if no effect beyond the 2a level is seen, then there will be a significant constraint on the 
model. 

The most recent experimental value of is averaged to be a e ^ p = 1165923(8.4) x 10~ 9 . 
The quantity within the parenthesis refers to the uncertainty in the last digit. The 
Standard Model results consist of several parts: 

a™ = < + < d + <~ W (1) 

Here a q ^ d |fj is computed to 0(a 5 ) QED corrections. a^ ad consists of 0(a 2 ) |i| [0| 
and 0(a 3 ) || hadronic vacuum polarizations, and also light-by-light hadronic contribu- 
tions ||21|| . We exhibit two different evaluations in Table 1. Case A uses the analysis of 
Martinovic and Dubnicka [|ll| who use a detailed structure of pion and kaon form factors 
in their fits to get the 0(a 2 ) corrections. Case B uses the 0(a 2 ) analysis of Eidelman and 



Jegerlehner fl9|, 120]. We note that there is almost a factor of two difference between the 



hadronic errors of case A and case B. For the electro-weak corrections we have used the 



recent analyses of Czarnecki et al. [22|, which include one- loop electro- weak corrections of 



4 



the Standard Model [P3"| , 2-loop corrections with fermion loops 23fl , and partial 2- loop 
corrections with boson loops [25] . The remaining unknown 2-loop corrections with boson 



loops are denoted by Rb following the notation of Ref. [p2 |. Pending the full 2-loop bosonic 
contributions in the E-W corrections, the total a^ M for case A is 11659180(7.6) xl0~ 10 
while for case B one has 11659200(15.3) xlO -10 . We see that the overall error, which is 
dominated by the hadronic corrections, is about half for case A relative to that for case B. 

The new high precision E821 Brookhaven experiment Q has an anticipated design 
sensitivity of 4 x 10~ 10 . This is about 20 times more accurate than that of the CERN 
measurement conducted earlier. However, as mentioned already one needs in addition an 
improvement in the computation of the hadronic contribution a l Jf ,d , where uncertainties 
primarily arise due to hadronic vacuum polarization effects. This problem is expected to 
be overcome soon through further accurate measurements of <T to4 (e + e~ — > hadrons) for 
the low energy domain. 

2 The Minimal Supergravity Model 



The framework of our analysis is N=l Supergravity grand unified theory W2\ where the 
Supergravity interactions spontaneously break supersymmetry at the Planck scale (Mpi = 
2.4 x 10 18 GeV) via a hidden sector fll2|. Further we assume that the GUT group G breaks 



at scale Q = M G to the Standard Model gauge group : G — ► SU(3) C x SU(2) L x U(1) Y - 
At low energy the symmetry breaking effective potential below the GUT scale Mq is given 
by 
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V SB = mlziz\ + (A W® + B W® + h.c.) (2) 

where W<® ,W<® are the quadratic and the cubic parts of the effective superpotential. 
Here W {2) = fx H 1 H 2 , with H u H 2 being two Higgs doublets, and contains terms 
cubic in fields and involves the interactions of quarks, leptons and Higgs with strength 
determined by Yukawa couplings. In addition the low energy theory has a universal 
gaugino mass term — mi\ a X a . The minimal Supergravity model below the GUT scale 
depends on the following set of parameters. 

m ,mi,A ,B ] fi ,a G ,M G (3) 

where M G and a G are the GUT mass and the coupling constant respectively. Among these 
parameters a G and M G are determined with the high precision LEP results by using two- 



loop renormalization group equations ai(Mz), i — 1, 2, 3 up to M G [26|. Renormalization 



group analysis is used to break the electro-weak symmetry |13[ and radiative breaking 

<fli>' 



allows one to determine //q by fixing Mz and to find Bq in terms of tanf3 = 7^4- Thus 



the model is completely parametrized by just four quantities [27 



m ,mg, A t , tan/3 (4) 
Here the universal gaugino mass mi is replaced by the gluino mass m„ and A by A t , 

2 a 

which is the t-quark A parameter at the electro- weak scale Mew- in addition to these 
four parameters and the top quark mass m f , one has to specify the sign of /i, since the 
radiative breaking equations determine only /i 2 . There are 32 new particles in this model 
(12 squarks, 9 sleptons, 2 charginos, 4 neutralinos, 1 gluino, 2 CP even neutral Higgs, 1 
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CP odd neutral Higgs and 1 charged Higgs). The masses of these 32 new particles and all 
their interactions can be determined by the four parameters mentioned above. Thus the 
theory makes many new predictions f27 |, and has led to considerable activity [ 28| , p9| to 
explore the implications of supergravity grand unification and its signals [EC]. The allowed 
parameter space of the model is further constrained by the following i) charge and colour 
conservation [^TJ, ii) absence of tachyonic particles, iii) a lower bound on SUSY particle 
masses as indicated by CDF, DO, and LEP data,iv) an upper limit on SUSY masses from 
the naturalness criterion which is taken as mo,mg < 1 TeV.Our analysis automatically 
takes into account important Landau pole effects that arise due the top being heavy 



and thus in proximity to the Landau pole that lies in the top Yukawa coupling |32 
Additionally, we consider the constraint on b — > 57 decay rate from the recent CLEO 
data, and the neutralino relic density constraint. 



3 Analysis of (g - 2)^ and Results 

We use Ref. || to compute SUSY contributions to (g — 2) M . These contributions arise from 
Figs, (la) and (lb). In Fig. (la) one exchanges 2 charginos W a , a=l,2 with masses m^ a 
which are charged spin ~ Dirac fields and a sneutrino state. In Fig. (lb) one exchanges 

4 neutralinos Z^, k=l,2,3,4 with masses m z {k) which are spin | Majorana fields (our 
labelling of particles satisfy m; < rhj for i < j), and two scalar smuon mass eigenstates. 
The mass spectra of the charginos, neutralinos, smuons and of the sneutrino is given in 
the Appendix for convenience. The one-loop supersymmetric contribution to (g — 2) M 
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then is given by 



,SUSY 



(5) 



For the chargino-sneutrino part referring to M we find (with summation over repeated 
indices implied) 



9* = ^k^-FiM + ^^^F 2 (x a ) 



247T 2 m\ 

W a 



4vr 2 



m 



(6) 



Wa 



a = 1,2 



m 



(7) 



W a . 



The mass of the sneutrino required in Eq. (pi) can be determined from Eq. fl33|) . The first 
term in Eq. (^) contains terms diagonal in chirality whereas the second term has R-L 
interference terms arising from Yukawa interactions. The functions F\(x) and F 2 (x) are 
defined by 



F x {x) = (l-5x-2x 2 )(l- x)~ 3 -6x 2 (l- x)~ 4 ln(x) 
F 2 (x) = (l-3x)(l-x)" 2 -2x 2 (l- x)~ 3 ln(x) 



A$ and A ( ^' of Eq. are given as 



| (a) 



.4 



(i) 

R 



(1) 



\/2 sin W 



cos 7i ; A\ 



m M cos 72 



a/2 sin W v^Mjy cos ,5 



(8) 



(9) 



.4 



(2) 
R 



sin 7a; 



(2) 



m„ sin 72 



a/2 sin 6^ y/2 sin 0yp \/2Mw cos (3 

Here the angles 71 and 72 are found by using 7 ij2 = /3 2 T Pi where 



sin 2/?i j2 
cos 2/3i 2 



(/i ± m 2 ) 4z/ 2 2 + (/i ± m 2 ) 



2^1,2 4i/ 2 2 + (// ± m 2 ) 



(10) 
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where — tt < 2/3i ;2 < 7r. 

The neutralino-smuon loop correction results in 
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K 1 

12tt 2 m| 



i(fc) 



5,X L [G 2 (a: lfc )-G 2 (x 2fc )] 



47r 2 



47r 2 



\0 k J2-dR 



{c*B L k - {-iy«s*B«}G,{x lk ) + {s 2 Bt - (-l) e «c 2 B«}G 2 (x 2k ) (12) 



Here s = sin 5, c = cos 5, and 



m 



x rk 



rrir 



r=l,2; fc = l,2,3,4 



4 (*o . 



The functions Gi(:r) and G 2 (a;) are given by 
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Gi(ar) = -(2 + 5x-a; 2 )(l-a;)~ 3 + 3a;(l-a;)" 4 /n(x) 
G 2 (x) = (l + x)(l-x)~ 2 + 2x(l-x)~ 3 /n(x) 



(13) 



(14) 



The Yukawa coefficients Cfc are found from 



Cfc = e- 



2Mi4/cos /3sin ^14/ 
and and B k are found by using 



[cos (30 3k + sin (30 Ak ] 



(15) 



5^ = e [-0 lk + cot 2^0 2fc ] 

5 fe L = e [-0 lk -t a n6 w 2k ](-l) 9k (16) 



The angle 5, 9 k and the quantities are defined in the Appendix 



Among the two sources of one-loop contributions to a^ USY the chargino-sneutrino 
loop contributes more than the neutralino-smuon loop. This occurs mainly due to the 
smallness of the mixing angle S (see Eq. (|12|)) which itself arises from the smallness of the 
muon to the sparticle mass ratio (see Eq. (|5TD). Partial cancellations on the right hand 
side of Eq. (O) are also responsible for a reduction of the neutralino-smuon contribution. 
For the chargino-sneutrino contribution (see Eq. @) one has comparable magnitudes for 
chirality diagonal and non-diagonal terms for small tan/3(~ 1). The non-diagonal terms 
become more important as tan/5 starts to deviate significantly from unity. Indeed, a large 
value of tan/3 results in a large contribution of the chirality non-diagonal term in the 
chargino-sneutrino part and hence to |a'? f/ ' sy | due to the enhancement arising from 

of I I o COS/ g 

~ tan/3 in the Yukawa coupling. 

There exists a very strong correlation between the sign of a^ c/sy and the sign of fi 
which we now explain. While the chiral interference chargino part of dominates 
over the other terms, it is the lighter chargino mass which contributes most dominantly. 
In fact this part depends on A^A^ which from Eq. (^) can be seen to have a front 
factor of (— l) e , where 6 = 0(1) for Ai > 0(< 0) where Ai is the smaller eigenvalue of the 
chargino mass matrix (see Eq. (|26|) ). For ji > 0, one finds Ai < invariably and for \i < 
one has Ai > for almost all the regions of parameter space of interest. This can be seen 
by writing Ai 2 in the following form 



1 



Ai, 2 = - (|2M^ + (/i - m 2 y - 2Mj v sin2(3\ 2 =f [2M^ + + m 2 y + 2Mf v sin2(3 

(17) 

and noting that the terms containing sin2/? are only appreciable when tan/3 is small (~ 1). 
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As a result, due to this unique dominance of the chiral interference term involving the 
lighter chargino mass, one finds as a rule a^ SY > for /i > and a^ USY < for fi < 
with some very few exceptions for the latter case when tan/3 ~ 1 resulting in a very small 

\ n SUSY\ 
\ a ti I- 

The generic dependence of |af c/s:r | on the remaining parameters m , m g , and A t is 
as follows: Regarding mo, the dependence results primarily from the chargino-sneutrino 
sector because the mass spectra depend on m$. This results in a decreasing |a^' 75y | with 
increase in mo. Among the other two basic parameters a large gluino mass m g in general 
again leads to a smaller |a^ c/5y |, due to the resulting larger sparticle masses entering in the 
loop. The A t dependence enters implicitly via the SUSY mass spectra, and also explicitly 
in the neutralino-smuon exchange diagrams. Figs. (2a) and (2b) show the upper limit of 
l a f USY \ in the minimal Supergravity model by mapping the entire parameter space for 
m o, m § < 1 TeV, tan/3 < 30 and both signs of fi. One finds that max \a^ USY \ increases 
with increasing tan/3 for fixed m g . As discussed above this happens due to the essentially 
linear dependence of the chirality non-diagonal term on tan/3. Furthermore, for large tan/? 
one finds similar magnitudes of max \a^ USY | for fi > and fi < 0. This can be understood 
by noting that for a large tan/3 the lighter chargino masses for /i > and < cases 
have almost similar magnitudes (see Eq. (|I7D). This, along with the explanation of the 
dependence of the sign of on the sign of fi accounts for similar la^ 5 ^ values for 

both fi < and fi > when tan/3 is large. This similarity between the fi > and fi < 
cases is less apparent for small tan/3 (ie. tan/3 ^ 2) and small m g values. 

Next we include in the analysis b — > 57 constraint and the dark matter constraint 
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which have been shown in recent work |33|, [34] to generate strong constraints on the 
parameter space. For the b — > 57 decay the CLEO Collaboration [17] gives a value of 



BR{b -> 57) = (2.32 ± 0.5 ± 0.29 ± 0.32) x 10" 4 (18) 
Combining the errors in quadrature one has BR(5 — > 57) = (2.32 ± 0.66) x 10~ 4 . 



The Standard Model prediction for this branching ratio has an 0(30%) uncertainty [35 
mainly from the currently unknown next-to-leading (NLO) order QCD corrections. Recent 
Standard Model evaluations give || BR(6 -> 37) = (2.9 ± 0.8) x 10~ 4 at m t » UOGeV. 
The SUSY effects in BR(6 — > sj) can be conveniently parametrized by introducing the 
parameter r S usY which we define by the ratio 



rsusr = BR(b -> s^) SUSY /BR{b -> sj) S m ■ (19) 

Several uncertainties that are present in the individual branching ratios cancel out in 
the ratio tsusy- However, we point out that the NLO corrections would in general be 
different for the SUSY case than for the SM case due to the presence of different SUSY 
thresholds []38] . In the present analysis we limit ourselves to the leading order evaluation. 
In a similar fashion we can define 

r exp = BR(b -> S1 ) exp /BR(b -> S1 ) SM • (20) 

Using the experimental result of Eq. (|T8j) and the Standard Model values given above one 
finds that r exp lies in the range 

r exp = 0.46 - 2.2 (21) 
12 



Now normally in SUSY theory one can get rather large deviations from the SM results so 
that tsusy can lie in a rather large range, i.e., ~ (0, 10). Thus the constraint tsusy = r exp 
is an important constraint on the theory. This constraint then excludes a part of the 
parameter space B3, |36| and reduces the magnitude of the max \a^ USY \ for a given m- g . 



The cosmological constraint on the neutralino relic density |5S] also plays a very sig- 
nificant role in limiting the parameter space of the model. Theoretically the quantity of 
interest is fi^/? 2 where Vt% i = -^ x . Here p^ x is the neutralino mass density, p c is the crit- 
ical mass density to close the universe and h = H/ (100km/ sMpc) where H is the Hubble 
constant. Astronomical observations indicate h = 0.4 — 0.8 which results in a spread of 
value for f2^ /i 2 . For our analysis we use a mixture of cold and hot dark matter (CHDM) 
in the ratio Qcdm '■ ^hx>a/=2:1 consistent with the COBE data. Assuming total Q = 1 
as is implied by the inflationary scenario and using for the baryonic matter f2#=0.1 we 



ret p| 3S 



0.1 <VLzh 2 <0A (22) 

The combined effects of the b — * 57 constraint and relic constraint put severe limits 
on the parameter space [36|, [ID, Hfl . Their effect on a^ USY is shown in Figs. (2c) and (2d) 



which are similar to Figs. (2a) and (2b) except for the inclusion of the combined effects of 
b — > 57 constraint and relic constraint. Comparison of Figs. (2a) and (2c) and similarly 
of Figs. (2b) and (2d) show that typically max |af | for large tan/3 ( i.e., tan/3 ^ 2) is 
reduced by about a factor of | for gluino masses below the dip under the combined effects 
of b — >• 57 and dark matter constraints. However, as is obvious, the most striking effect 
arises due to the appearance of the dip itself. The existence of such a dip was first noted in 
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Ref. |36| and is due to the relic density constraint. It is caused by the rapid annihilation 
of neutralinos near the Z pole which reduces the relic density below the lower limit in 
Eq. (|22|) and hence part of the parameter space gets eliminated due to this constraint. In 
order to get the correct position and depth for this dip one must use the accurate method 
3T| for calculation of the relic density. The analysis of Figs. (2c) and (2d) show that 



it would be difficult to discern SUSY effects in the (g — 2) M experiment for gluino masses 
around the dip or correspondingly for neutralino masses of m^ ~ Mz/2. ( There is a 
similar dip at ~ m/j/2 which does not appear in the graphs because the remaining 
parameters have been allowed to range over the full space , but would be manifest once 
the Higgs mass is fixed.) 

Interestingly, even the current limits on (g — 2) M including the present experimental 
and theoretical errors place some constraint on the parameter space of supergravity grand 
unification. Ascribing any new physics to supersymmetry by using + a ^ M — a ^T P one 

may constrain the parameter space of the model in the (mo — Trig) plane for different tan/3 
with a consideration of all possible A t . We have combined the uncertainty of theoretical 
estimate and the current experimental uncertainties in quadrature to find that 



- 11.7 x 1(T 9 < al USY < 22.5 x 1CT 9 (23) 

Fig. (2e) exhibits for the case < the excluded regions in the (mo — trig) plane 
for different tan/3 values where the excluded domains lie below the curves. We note that 
the excluded domains depend strongly on the value of tan/3 and constraints become more 
severe as tan/3 increases. For tan/3 ^ 10 the constraints on (mo,mp) are very modest. A 
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similar analysis holds for /i > as shown in Fig. (2f), except that here the constraints on 
(mo, rrig) are generally less severe. 

We analyze now the effect of the implications of the predicted experimental accuracy 
of (~ 4 x 10 -10 ) to be attained in the Brookhaven experiment for supergravity grand 
unification. Of course, the present theoretical uncertainty in the hadronic contributions 
to a M , mostly arising from the hadronic vacuum polarization effects as discussed earlier, 
limits the usefulness of such a precise measurement of a^. The hadronic uncertainty arises 
from the uncertainty in the computation of [|J 

D , x a to i(e + e~ -> hadrons) 

H \ s ) = rrr^ T~ \ — I 24 ) 

cr(e+e — > u+fi ) 

for the low energy tail of (e + e~ — > hadrons) cross-section. Ongoing measurements in 
VEPP-2M together with future experiments in DA$NE, BEPC etc. are expected to 
reduce this hadronic uncertainty to a significant extent, perhaps by a factor of 4 or more, 
enhancing correspondingly the usefulness of the precision measurement of a^. 

We present here two analyses. For the first analysis we used the more optimistic 
estimate of Ref. [IIS) where the authors made an improved evaluation of R(s) and a^ ad 



through the use of global analytical models of pion and kaon form factors in addition 
to the use of a better experimental input of the three-pion e + e~ annihilation data in 
comparison to previous determinations [[J. This corresponds to Case A of Table 1 and 
gives (setting Rb = 0) the result 

af 1 = 11659183(7.6) x 1(T 10 (25) 

For the second analysis we shall make a comparative study over an assumed range of 
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errors which includes analyses of both cases A and B of Table 1 as well as the possibility 
of even more constrained errors. 

In order to analyze the effect of the predicted accuracy level of the Brookhaven ex- 
periment on the models of our discussion we have assumed that the experiment will not 
discard the Standard Model result within its 2a uncertainty limit. As in the analysis 
of the current experiment above we ascribe any new physics to supersymmetry by using 
a susY _^ a SM _ a exp an( j cons t ra j n the parameter space of our model in the (mo— m^) plane 
for different tan/3 with a consideration of all possible A t . Following the same procedure as 
before, we combine the uncertainty of theoretical estimate and the expected experimental 
accuracy level of 4 x 10~ 10 in quadrature. In our first analysis we use Eq. ( ^5|) as the 
theoretical input and assume that the predicted accuracy in the experimental determi- 
nation of a M will be achieved, we then determine the constraints on the SUSY particle 
spectrum if a^ USY lies within the 2a of the combined theoretical and experimental error. 
In Figs. (3a) to (3f) we exhibit the regions of (m^ — Trip ) plane which will be excluded 
(dark shaded region) if a^ USY lies within the 2a limit. We also exhibit the regions which 
will be partially excluded (light shaded area) because a significant part of the parameter 
space is eliminated by the constraint, and the allowed region (dotted area). In addition, 
there is a region (white space) which is inaccessible due to the existence of a lower limit 
of sneutrino mass (m? ) [j36fl . 

Next we give a comparative analysis of the constraints for cases A and B of Table 1 
and also for a case C where the error is reduced by factor of 4 over case B ( see in this 



context the analysis of Ref. [57] which gives a new evaluation of hadronic contributions 



16 



of 699(4.5) x 10~ 10 ). In Figs. (4a) and (4b), we give a composite display of the excluded 
regions in the (mo — rrig) plane for the three cases. We observe that the forbidden region 
increases in proportion to the decrease in the combined error of theory and experiment. 
Further, the excluded region increases with increasing value of tan/3. Thus if the combined 
error decreases by a factor of four (case C) and no effect beyond la is seen, then the g — 2 
experiment will exclude most of the region in (mo — rrig) plane for large tan/3 as can be 
seen from Figs. (4e) and (4f). In fact even with the presently large uncertainty in the 
theoretical values of a^ M one will be able to exclude a significant part of the (mo — rrig) 
parameter space if the expected sensitivity of the measurement of a M is reached for large 
tan/3 values (see Figs. (4e) and (4f)). Of course a significant reduction of the uncertainty 
in a^ M which one expects to be possible in the near future will more stringently constrain 
the parameter space when combined with the expected sensitivity of the Brookhaven 
experiment. We have also carried out a similar analysis for the af USY for the tau lepton. 
This gives max \a^ USY \ ~ 1.0 x 10~ 5 for tan/3 = 20. Even for this large value of tan/3 the 
predicted value of af USY is too small to observe with present experimental accuracy |43| . 



4 Conclusion 

In this paper we have presented an analysis of (g — 2)^ USY within the framework of 
supergravity grand unification under the constraint of radiative breaking of the electro- 
weak symmetry, and the constraint of b — > S7 and dark matter. One finds that over most 
of the parameter space the chiral interference light chargino part of a^ USY dominates 
and imparts a signature to a^ USY .Thus as a rule one finds a^ USY > for ii > and 
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a jF < for < with very few exceptions when tan/5 ~ 1. At the quantitative 
level it is shown that the expected experimental sensitivity of (g — 2) M measurements 
combined with the expected reduction of error in (g — 2)^ M ( by pa 0(1/4)) will exclude 
the m ,mg parameters in the domain m < 700 GeV, rrig < 1 TeV for large tan/3(> 20) 
and stringently constrain the parameter space for lower tan/3. With the same assumptions 
one will be able to probe the domain m ^ 400 GeV, m- g ^ 700 GeV for tan/5 ^ 10. The 
constraint becomes less severe for smaller values of tan/3. However, even for tan/3 = 5, the 
SUSY reach of the new experiment will be very substantial (see Figs. (4a) and (4b)). Thus 
one finds that the Brookhaven experiment coupled with the corroborating experiments 
and analyses designed to reduce the hadronic error will complement SUSY searches at 
colliders and provide an important probe of the parameter space of supergravity grand 
unification especially for large tan/3. 

We wish to acknowledge useful discussions with R. Arnowitt, W. Marciano and W. Worstell. 
This research was supported in part by NSF grant number PHY-19306906 and PHY94- 
07194. 



5 Appendix 

The chargino masses m\y = |Aj|, i=l ,2 where \\^ are the eigenvalues of the chargino 
mass matrix and are given by 



A l,2 = \{ ~ (/' - '''2)" 



Av\ + (/i + m 2 y 



(26) 
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where 



Mi 



^1,2 



72 



sin (3 =F cos /3) 



(27) 



and 777.2 is obtained from the relation 777 a 



a a (M z ) 
a z {M z ) 



nig where a a (Mz), a = 1,2,3 are the 



SU(3), SU(2) and U(l) gauge coupling constants at the Z boson mass. 



The neutralino masses m 



Z(k) 



|Afc| where Xk are the eigenvalues of the neutralino 



mass matrix which in the (W3, B, Hi, H 2 ) basis reads [12 



M z 



1 m 2 





a 


b 





fhi 


c 


d 


a 


c 







V b 


d 


-n 






(2f 



where a = Mz cos 6w cos j3, b 
Mz sin 9w sin (3 while the quantities 6k that appear in Eq. fll6|) are defined by 6k = (1) 
for Ajt > (< 0). The quantities are the elements of the orthogonal matrix which 
diagonalizes the neutralino mass matrix. 
Smuon masses are given by 



777 ,- 



-{(ml + rh 2 R ) =F sj {m\ - m\) 2 + Am^ (A t + // cot ^ 



M1.2 2 ' 

Here is scaled with mo, and the L and R parts are given by 



(29) 



777? 



777p 



ml + m 2 + a G 



3 f + 3 f 



ml + + sin z 6 W ) M z z cos 2/3 



m o + m M + ( ^ /il m ± — s i n2 @wMz cos 2/3 



(30) 



where mi = ^^ym^, a G = a G /4vr, / fc (*) = t(2+/3 fc 7;)/(l+/? fe t) 2 with/3 fc = a G (33/5, 1, -3) 
and t = 21n(^p) at Q = M z . 
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The mixing angle which describes the left-right mixing for the smuons is determined 
by the relations 

• ox 2m lx (A t + /icot/3) 

sin 2d = , (31) 
^{ml-rhlf + Aml (A t + fj, cot (3) 2 



and 



2m M (A t + /i cot (3) 

tan 2d = i 7-n 5 (32) 

{ml - m\) 



A similar result analysis holds for sneutrino masses and one has 



m? M = ml + a G 



3 3 
2 /2 + 10 71 



mi 2 + iM|cos2/5 (33) 
2 2 
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Figure Captions 

Fig. 1: 

Fig. (la) Chargino-sneutrino one-loop exchange diagram which contributes to o£ . 
Fig. (lb) Neutralino-smuons one- loop exchange diagram which contributes to a^ SY . 



Fig. 2: 

Fig. (2a) The upper limit of |af | vs. m g for the case \i < in the Minimal Supergravity 
model for different values of tan/3 in the region tan/3 < 30, when one allows the remain- 
ing parameters A t and mo to vary over the parameter space subject to the naturalness 
criterion of m < 1 TeV. 

The dashed horizontal line is the current 2a limit as given by Eq. (p3|) 
Fig. (2b): Same as Fig. (2a) except that u, > 0. 

Fig. (2c): Same as Fig. (2a) except that the constraint from b — > sj and dark matter as 

discussed in the text are included. 

Fig. (2d): Same as Fig. (2c) except that // > 0. 

Fig. (2e): Excluded regions (regions enclosed by the curves towards the origin) in the 
(m — nig) plane under the current theoretical and experimental limits of a M for various 
values of tan/3, for the case \i < . 
Fig. (2f): Same as Fig. (2e) except that /i > 0. 



Fig. 3: 

Fig. (3a) Display in the (ttw — m^) plane of the allowed (shown in dots), disallowed 
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(shown in squares) and partially allowed (shown in cross) regions corresponding to the 2a 

limit of Case A of the theoretical evaluation of a M and the predicted level of accuracy of 

the Brookhaven experiment, for the case tan/5 = 5 and /i < 0. The white area between 

the excluded region and axis remains inaccessible due to the existence of a lower 

limit of sneutrino mass rrip . 

Fig. (3b): Same as Fig. (3a) except that fi > 0. 

Fig. (3c): Same as Fig. (3a) except that tan/? = 10. 

Fig. (3d): Same as Fig. (3c) except that \i > 0. 

Fig. (3e): Same as Fig. (3c) except that tan/5 = 30. 

Fig. (3f): Same as Fig. (3e) except that /i > 0. 



Fig. 4: 

Fig. (4a): Excluded regions (regions enclosed by the curves towards the origin) for 
tan/? = 5 and < in the 2a limit, after combining in quadrature the predicted 
Brookhaven experimental uncertainty and different levels of uncertainty of theoretical 
estimates corresponding to Cases A, B, and C as discussed in the text. 
Fig. (4b): Same as Fig. (4a) except that ji > 0. 
Fig. (4c): Same as Fig. (4a) except that tan/5=10. 
Fig. (4d): Same as Fig. (4c) except that /i > 0. 
Fig. (4e): Same as Fig. (4a) except that tan/3=30. 
Fig. (4f): Same as Fig. (4e) except that /i > 0. 
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Fig. (la) 




Fig. (lb) 
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